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$8 3e. 根据 土壤 盐 涡 化 空间 异 质 性 对 南 疆 干旱 区 绿洲 农田 进行 精准 管理 分 区 ,对 农业 种 植 结构 调整 和 精细 化 管 


理 具 有 重要 意义 。 以 典型 干旱 区 绿洲 农田 为 研究 对 象 ,以 电磁 感应 数据 .地 形 数据 和 卫星 遥感 数据 为 数据 源 ,通过 


地 统计 学 方法 分 析 研 究 区 土壤 盐 渍 化 的 空间 异 质 性 ,利用 相关 分 析 筛 选 出 不 同时 


期 的 植被 指数 和 盐分 指数 。 以 农 


田 表 层 表 观 电导 率 (ECiosas) 为 主 变量 ,深层 表 观 电导 率 (ECos ECoxs ECs) .植被 指数 (RVI、GRVI、EVI) .土壤 盐分 


指数 (NDSI.S5 .SILT ) 为 辅助 变量 ,采用 
Coefficient of Variation ) 和 莫 兰 指数 


面向 对 象 多 尺度 分 割 算法 对 研究 区 进行 分 区 ,并 用 分 区 内 平均 变异 系数 (Cv， 
Moran's 1) 对 分 区 结果 进行 评价 。 结 果 表 明 :(1) 研究 区 各 土 层 表 观 电导 率 都 


存在 明显 空间 异 质 性 , 且 各 辅助 变量 与 主 变 量 均 存在 极 显 著 相关 。(2) 各 分 区 的 平均 变异 系数 较 整个 研究 区 的 变异 


系数 下 降 了 60%, 且 基于 多 源 数据 分 区 的 区 间 异 质 性 均 高 于 单一 数据 分 区 的 区 间 异 质 性 ,(3) 从 农田 耕作 
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农田 管理 又 符合 机 械 化 作业 要 求 。 人 研究 结果 可 为 南 疆 干旱 区 绿洲 农田 的 精准 管理 分 区 提供 一 定 的 技术 和 方法 


借鉴 。 


关键 词 : 土壤 盐 渍 化 ; 多 源 数据 ; 表 观 电导 率 ; 面向 对 象 分 割 ; 精准 管理 分 区 ; 


土壤 盐 渍 化 不 仅 是 全 球面 临 的 生态 问题 " , 同 
时 也 是 影响 全 球 灌溉 农业 最 重要 的 因素 之 一 ,而 
且 不 同 程度 盐 溃 化 土壤 对 不 同 植物 和 农作物 生长 
发 育 影响 也 不 尽 相 同 。 因 此 ,掌握 土壤 盐 省 化 信息 
可 在 一 定 程度 上 指导 农业 合理 灌溉 和 合理 调整 作 
物 布局 。 近 年 来 , 随 气候 的 变化 ,特别 是 全 球 气候 
变 暖 将 显著 影响 全 球 水 循环 ,导致 农业 用 水 量 急剧 
减少 ”。 在 干旱 区 . 半 干 旱 区 ,土壤 盐 渍 化 问题 和 农 
业 汐 溉 缺 水 问题 尤为 突出 。 因 此 在 干旱 区 EF 
量 区 急需 一 种 根据 盐 渍 化 程度 进行 变量 灌溉 和 调 
整 作 物 布局 的 技术 ,而 精准 管理 分 区 是 解决 此 类 
问题 的 根本 途径 ,同时 也 是 当前 精准 农业 研究 的 
热点 之 一 “。 
精准 管理 分 区 是 根据 土壤 中 各 种 制约 产量 的 
影响 因子 ,划分 出 具有 相同 均 质 的 子 区 域 ,来 调整 
农业 投入 品 和 作物 布局 的 措施 ”“。 分 区 之 间 实 施 
差异 化 管理 措施 有 利于 提高 农业 生产 投入 效率 、 市 


收 稿 日 期 : 2021-08-21; ”修订 日 期 : 2021-12-03 


干旱 区 


约 生产 成 本 保护 生态 环境 ,从 而 实现 农业 的 可 持 
续 发 展 。 目 前 ,在 精准 管理 分 区 方面 ,国内 外 学 者 
大 多 研究 采用 土壤 养分 数据 和 产量 相关 数据 进行 
插值 ,再 利用 模糊 聚 类 算法 或 自然 断 点 法 进行 分 区 
管理 ,但 该 方法 在 数据 获取 方面 需要 耗费 大 量 的 
人 力 和 物力 , 且 不 易 大 面积 推广 。 近 年 来 , 随 着 于 
感 技 术 的 发 展 ,许多 学 者 开始 借助 遥感 影像 数据 ， 
再 利用 面向 对 象 多 太 度 分 割 算法 或 模糊 聚 类 等 算 
法 进行 分 区 管理 "”, 但 此 方法 只 能 间接 的 通过 光 
谱 指数 反映 出 土壤 表层 属性 数据 ,而 不 能 直接 获取 
到 深层 土壤 属性 数据 。 同 时 也 有 部 分 学 者 使 用 土 
二 属性 数据 与 遥感 影像 数据 相 结 合 的 多 源 数据 , 利 
用 面向 对 象 多 尺度 分 割 算法 进行 分 区 管理 ,但 该 
方法 在 土壤 属性 数据 的 获取 方面 也 存在 费时 又 费 
力 的 缺点 ,只 能 在 较 小 区 域 展开 研究 。 上 述 研究 均 
是 针对 土壤 养分 差异 进行 管理 分 区 ,基于 土壤 盐 注 
化 差异 进行 管理 分 区 的 研究 其 少 。 干旱 区 或 半 干 
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旱 区 ,特别 是 新 疆 的 南部 地 区 ,土壤 盐 渍 化 是 影响 
作物 产量 和 品质 的 关键 因子 , 需 根 据 土 壤 盐 渍 化 程 
度 采取 不 同 的 管理 措施 ,如 作物 种 植 结构 的 合理 布 
局 . 冬 灌 或 春灌 的 灌溉 量 定额 。 因 此 在 这 种 土壤 盐 
渍 化 问题 突出 且 盐 渍 化 变异 强 的 南 疆 区 域 , 应 以 土 
壤 盐 渍 化 差异 作为 分 区 关键 因子 进行 精准 管理 分 
区 ,才能 提高 农业 效益 。 传 统 的 土壤 盐 渍 化 信息 获 
取 方 法 主要 依靠 田间 调查 结合 室内 分 析 ,存在 费 
时 、 费 力 ` 成 本 高 等 方面 的 不 足 , 再 考虑 到 盐 渍 化 强 
烈 的 时 空 变异 特征 ,该 方法 很 难 满足 区 域 尺 度 土壤 
盐 渍 化 的 调查 ,而 以 电磁 感应 为 原理 的 EM38-MK2 
大 地 电导 率 仪 ,可 快速 .高效 地 直接 获取 不 同 深度 
土壤 的 表 观 电导 率 数据 ,土壤 表 观 电导 率 数据 与 土 
壤 盐 分 含量 具有 高 度 正 相 关 性 ,可 有 效 反 映 土壤 
盐 渍 化 程度 ,使 用 EM38-MK2 大 地 电导 率 仪 可 显著 
提高 土壤 盐 渍 化 信息 的 获取 效率 ,特别 是 针对 区 域 
尺度 的 调查 而 言 。 自 20 世纪 末 以 来 ,使 用 表 观 电导 
率 来 表征 土壤 盐 涡 化 的 方法 在 土壤 盐 渍 化 监测 和 
评价 中 已 得 到 了 广泛 应 用 "” ,但 直接 应 用 表 观 电 
导 率 数据 进行 农田 管理 分 区 的 报道 鲜 见 。 此 外 , 目 
前 关于 农田 管理 分 区 的 研究 普遍 只 采用 了 土壤 的 
基本 信息 ,缺乏 综合 利用 土壤 信息 和 作物 信息 进行 
管理 分 区 的 报道 ,而 卫星 遥感 技术 可 以 快速 获取 大 
面积 作物 的 动态 信息 ,为 这 种 研究 思路 的 实施 提供 
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了 有 力 条 件 。 因 此 可 使 用 土壤 表 观 电导 率 和 遥感 
影像 数据 进行 土壤 盐 涡 化 的 分 区 管理 。 

综 上 所 述 ,本 研究 以 南 疆 干旱 区 典型 绿洲 阿拉 
尔 必 区 农田 为 研究 对 象 , 以 研究 区 不 同 深 度 土壤 的 
表 观 电导 率 数 据 、 地 形 数 据 和 卫星 还 感 影像 数据 为 
数据 源 ,以 土壤 表层 的 表 观 电导 率 数据 为 主 变量 ， 
在 分 析 土 壤 表 观 电 导 率 的 空间 异 质 性 和 空间 自 相 
关 性 的 基础 上 ,通过 相关 分 析 筛 选 出 与 表 观 电导 达 
到 极 显 著 相 关 的 卫星 遥感 光谱 指数 ,并 将 其 作为 畏 
助 变量 ,再 利用 面向 对 象 多 尺度 分 割 算法 进行 不 同 
输入 量 下 土壤 盐 渍 化 的 分 区 ,并 对 不 同 输入 量 下 的 
管理 分 区 结果 进行 评价 对 比 ,找到 最 优 管理 分 区 , 
并 提取 出 研究 区 各 田 块 的 边界 ,从 而 根据 田 块 所 在 
的 盐 渍 化 分 区 指导 农田 的 精准 灌溉 和 作物 的 合理 
布局 ,以 实现 作物 优质 高 产 的 目的 。 


1 研究 区 概况 


研究 区 位 于 新 疆 维吾尔 族 自 治 区 一 师 阿拉 
尔 市 十 二 团 境内 (图 1) ,地 理 坐 标 为 81"19'31" 工 ， 
40?29'52"N ,面积 为 1144 hm ,地 处 塔里木 盆地 西北 
部 、 塔 克拉 玛 干 沙漠 北 缘 , 天 山南 蕴 。 境 内 主要 有 
叶 尔 羌 河 、 和 田 河 、 阿 克 苏 河和 塔里木 河 四 大 水 
Ro 冬季 严 寒 , 夏 季 炎 热 ,日 照 时 间 长 ,降水 稀少 ， 
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图 1 研究 区 位 置 及 样 点 分 布 示意 图 


Fig. 1 Location ofthe study area and distribution of sample points 
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蒸发 强烈 。 该 区 多 年 平均 降水 量 为 48.5 mm, 多 年 
平均 蒸发 量 为 1988 mm, 28 KE LEIA 40.99 , 3 HL AY 
Y IB i BR rg XC Fi EET A re Ux o DEAE DO SA 
平坦 ,土壤 质地 以 砂 坏 土 为 主 ,区 域内 农田 主要 种 
植 棉花 ,棉花 生育 期 内 均 采 用 膜 下 滴灌 方式 进行 灌 
溉 ,而 春灌 和 冬 灌 则 采用 大 水 漫灌 的 方式 。 大 水 漫 
灌 旨 在 将 整个 棉花 生育 期 内 运 移 到 土壤 表层 的 盐 
分 淋 洗 到 耕 层 以 下 ,以 此 降低 下 茬 作物 种 植 的 盐 害 
胁迫 ,其 灌溉 用 水 主要 来 自 塔里木 河 。 研 究 区 土壤 
盐 渍 化 程度 差异 较 大 ,种植 结构 单一 , 冬 春灌 采用 
统一 灌溉 量 进 行 灌溉 ,不仅 易 造 成 局 部 田 块 压 盐 不 
彻底 导致 盐 害 现象 产生 或 过 度 淋 洗 浪费 水 资源 , 同 
时 也 造成 研究 区 作物 布局 不 合理 导致 作物 经 济 效 
益 下 降 , 因 此 ,研究 区 迫切 需要 根据 不 同 盐 渍 化 程 
度 来 进行 变量 灌溉 和 种 植 结构 的 调整 。 


2 数据 的 获取 与 研究 方法 


2.1 多 源 数据 的 获取 

2.1.1 电磁 感应 数据 的 获取 ”土壤 表 观 电导 率 数据 
SA BUR CK Hl HL Sj [XL (EM38- MK2) ,在 棉花 收 
获 后 冬 灌 前 的 11 月 上 旬 进 行 土壤 表 观 电导 率 数据 
的 采集 。 每 次 采集 表 观 电导 率 时 ,首先 将 大 地 电导 
率 (EM38-MK2) 设 置 为 手动 测定 模式 ,并 对 仪器 进 
行 预 热 和 归 零 校 验 , 竺 仪器 校 验 完毕 后 ,进行 表 观 
电导 率 的 获取 ;在 每 个 样 点 采集 时 ,仪器 将 自动 记 
录 每 个 样 点 的 地 理 坐 标 ,每 个 采样 点 同时 采集 垂直 
模式 和 水 平 模式 两 套数 据 ;水平 模 式 采集 的 为 0~ 
0.375 m 和 0~0.75 m 的 2 个 深度 土 层 的 表 观 电导 率 
数据 , 即 ECoas 和 了 ECox, 垂 直 模 式 则 为 0~0.75 m 和 
0~1.5 了 的 两 个 深度 土 层 表 观 电导 率 数 据 , 即 下 Cr 
和 Cs, 相 邻 采 样 点 间距 小 于 300 m。 研 究 区 共 随 
机 均匀 采集 约 1200 组 表 观 电导 率 数 据 ( 图 1) ,采集 
数据 描述 性 统计 见 表 1。 

2.2.2 还 感 影像 数据 和 地 形 数据 的 获取 ”研究 采用 
来 源 于 美国 地 质 勘 查 局 (USGS, http://glovis.usgs. 


gov) 的 Landsat8 遥感 影像 ,卫星 重 访 周 期 为 16 d, 空 
间 分 辩 率 为 30 m。 本 文选 取 2020 年 棉花 生长 旺盛 
期 (9 月 5 日 ) 的 影像 和 棉花 收获 冬 灌 前 (11 月 4 日 ) 
的 裸 土 影像 。 影 像 经 辐射 定 标 、 大 气 校正 图像 裁 
前 等 预 处 理 后 ,使 用 ENVI 5.3 软件 计算 得 出 9 月 5 
日 影像 的 9 个 植被 指数 (RVI、DVI、GDVI、GRVI、ND- 
VIEVINNIR NRNG) 和 11 月 4 日 影像 的 13 个 土 
壤 盐 分 指数 (NDSI、S1、S2、S3、S5、S6、SI-T、SI.、SIl、 
SI2, SI3, Intl, Int2) 。 同 时 通过 地 理 空间 数据 云 
(http://www.gscloud.cn) F R 25 [i] 27 HEX Jy 30 m 的 
DEM Zitfz , fii JH ArcGIS 10.7 软 件 完 成 洼地 填充 . 影 
像 裁剪 .几何 精度 校正 和 高 程 数据 提取 。 

2.2 ”分割 方法 

22. 基于 面向 对 象 的 多 尺度 分 割 ”面向 对 象 多 尺 
度 分 割 是 影像 对 象 的 一 个 两 两 合并 的 过 程 ,合并 从 
均匀 分 布 在 影像 上 的 影像 对 象 开始 。 对 于 每 个 
影像 对 象 ,找到 一 个 使 它们 之 间 的 异 质 性 变化 最 小 
相 邻 影像 对 象 ,在 异 质 性 变化 最 小 处 ,2 个 对 象 合并 
成 1 个 更 大 的 对 象 ,每 个 对 象 在 每 个 循环 中 处 理 1 
次 ,直到 分 割 对 象 的 异 质 性 超过 用 户 给 定 的 阔 值 ， 
合并 过 程 停止 。 本 研究 尝试 基于 多 源 数据 的 土 
壤 盐 涡 化 空间 变异 进行 分 区 , TE eCognition 软件 中 
首先 确定 合适 的 形状 因子 (Shape) 和 紧 致 度 因子 
(Compactness ) ,然后 通过 调节 不 同 的 分 割 尺 度 参数 
(Scale parameter) ,使 分 区 之 间 表 观 电 导 率 异 质 性 最 
大 ,分 区 内 部 电导 率 同 质 性 最 大 。 

2.2.2 确定 最 优 管理 分 区 ”由 于 面向 对 象 分 割 的 原 
则 为 尽量 使 分 割 对 象 之 间 的 异 质 性 高 ,同时 使 分 割 
对 象 内 部 的 同 质 性 高 ,所 以 为 提高 分 区 内 部 的 同 
质 性 和 分 区 之 间 的 异 质 性 需要 寻找 最 优 分 区 尺度 ; 
根据 学 者 研究 ,可 使 用 平均 分 割 评价 指数 的 变化 
来 确定 最 优 分 制 尺 度 ,将 平均 分 割 评价 指数 和 尺 
度 参数 绘制 成 曲线 图 ,在 平均 分 割 评价 指数 开始 
产生 明显 变化 时 的 前 一 个 尺度 参数 为 最 优 分 割 尺 
度 参数 ""。 通 过 在 多 尺度 分 割 算法 中 调节 分 割 尺 


表 1 不 同 深度 表 观 电导 率 描述 


Tab.1 Descriptive statistics of apparent electrical conductivity at different depths 


表 观 电导 率 T] MÉU/(mS - m^?) 最 大 值 /(mS:m') 平均 值 /(mS.m') 变异 系数 /% 峰 度 偏 度 
ECioss 2.39 300.98 59.39 81.82 1.54 1.87 
ECioss 3.23 287.15 55.58 71.95 1.32 1.83 
EC.o; 1.17 280.04 54.67 74.52 1.16 1.72 
ECs 2.24 290.43 54.19 78.97 2.82 1.02 
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度 参数 将 研究 区 分 为 严 个 区 域 ,并 计算 出 不 同 分 割 
尺度 下 各 区 域 的 标准 差 和 各 区 域 与 邻 域 的 平均 差 
分 绝对 值 ,其 同 质 性 指数 用 标准 差 ( 6 ) 来 表示 , 5 
质 性 指数 用 平均 差分 绝对 值 ( AC, ) 来 表示 ,计算 如 
公式 (1)、(2) 所 示 ; 同 时 构建 对 象 的 分 割 评价 指数 
(SEI) 和 平均 分 割 评价 指标 (45E7) ,计算 如 公式 (3) 
和 (4) 所 示 。 


1 x 2 
ò= LL (6,- c) (1) 
1 N 
AC,» FEM (2) 
会 
SEI= Ber (3) 
8 
_ly 4. 
ASEI- 9 7 A; SEI, (4) 


式 中 ;nn 为 菜 区 域内 像 元 的 个 数 ; Ci 为 该 区 域内 第 
i 个 像 元 的 属性 值 ; C, 为 该 区 域 的 属性 平均 值 ; 6 表 
示 该 区 域 的 标准 差 ;了 工 为 该 区 域 的 周 长 ; N 表示 与 
该 区 域 相 邻 区 域 的 个 数 ; L 表示 与 第 ;个 相 邻 区 域 
公共 边 的 长 度 ; Cp 为 第 了】 个 相 邻 区 域 的 属性 平均 
值 ; AC, 表示 该 区 域 与 邻 域 的 平均 差分 绝对 值 ; SEI 
表示 分 割 评价 指数 ; A 表示 研究 区 的 总 面积 , 4, 为 i 
个 区 域 的 面积 ; m 为 被 划分 的 区 域 个 数 ; SET, 为 第 i 
个 区 域 的 分 割 评价 指数 ; ASEI 表示 平均 分 割 评价 
指数 。 

2.3 精准 管理 分 区 的 评价 方法 

2.3.1 分 区 内 部 的 同 质 性 ”根据 分 区 内 像 元 属性 值 
的 标准 差 和 平均 值 ,分 别 计算 出 各 个 分 区 内 的 变异 
系数 ,并 考虑 到 各 分 区 面积 大 小 不 同 而 赋予 相应 的 
权重 ,从 而 用 分 区 内 部 变异 系数 (Cv) 的 平均 值 表示 
同 质 性 ,计算 如 公式 (5) 所 示 : 


C, 2 =! 10096 (5) 


式 中 : Cy; 为 第 i 个 分 区 的 变异 系数 ; a; 为 第 ;个 分 
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区 的 面积 ; m 为 被 划分 的 区 域 个 数 。 C 越 小 ,分 区 
内 变异 越 小 ,代表 分 区 内 异 质 性 较 小 , 同 质 性 较 大 。 
2.3.2 分 区 之 间 的 异 质 性 采用 全 局 莫 兰 指数 ( 了) 
判断 空间 相关 性 ,根据 英 兰 指数 大 小 判断 分 区 间 蜡 
质 性 高 低 ,其 计算 如 公式 (6) 所 示 : 


DAC -¥\(y, -7) 
和 1 


Is 
(Bon Zr] 

式 中 : m 为 被 划分 的 区 域 个 数 ; y, 为 分 区 i 的 属性 
平均 值 ; 7 为 整个 研究 区 的 属性 平均 值 ; W, 为 i 和 j 
之 间 的 空间 权重 , 若 :与 了 相 邻 , 则 栈 =1 ,和 否则 
本 =0 。 了 表示 空间 相关 程度 , 7 越 接近 于 0, 表 示 
分 区 之 间 相 关 程 度 越 低 , 即 分 区 之 间 蜡 质 性 越 高 ， 
而 根据 分 区 评价 原则 可 知 ,分 区 之 间 异 质 性 越 高 ， 
分 区 结果 越 好 。 


3 结果 与 分 析 


3.1 土壤 盐 渍 化 的 空间 异 质 性 

根据 土壤 盐 渍 化 来 划分 管理 分 区 , 需 证 明 该 研 
究 区 土壤 表 观 电导 率 具 有 空间 异 质 性 ,因此 需 使 用 
地 统计 学 中 半 变 异 函 数 进行 分 析 。 在 计算 半 变 异 
函数 时 ,数据 需 符合 正 态 分 布 ,否则 会 出 现 比例 效 
应 。 经 分 析 , 原 始 数 据 呈 偏 态 分 布 ,经 平方 根 转化 
后 基本 符合 正 态 分 布 。 通 过 GS+ 软 件 获 得 变异 函数 
的 最 优 模型 ,其 参数 如 表 2 所 示 。 

由 表 2 可 知 ,各 个 土 层 表 观 电导 率 的 最 优 模型 
均 为 指数 模型 ,决定 系数 均 大 于 0.93 , 且 残 差 均 小 于 
1.27 ,说明 该 模型 可 较 好 地 反映 研究 区 土壤 盐 渍 化 
的 空间 分 布 情况 。 块 金 值 表示 人 为 随机 因素 造成 
的 空间 异 质 性 ,各 土 层 表 观 电导 率 的 块 金 值 在 3.05~ 
5.33 之 间 ,表明 研究 区 表 观 电导 率 存 在 由 人 为 随机 
因素 所 造成 的 空间 蜡 质 性 。 基 台 值 表示 人 为 随机 
因素 和 自然 结构 因素 造成 的 最 大 程度 空间 异 质 性 ， 


(6) 


表 2 不 同 深度 表 观 电导 率 空间 异 质 性 
Tab.2 Spatial heterogeneity analysis of apparent electrical conductivity at different depths 


表 观 电导 率 模型 块 金 值 基 台 值 块 金 值 / 基 台 值 /% 变 程 /m 决定 系数 残 差 
RC Ep) 5.33 10.79 49.38 2579.00 0.96 0.92 
ECros 指数 3.59 7.18 49.99 2148.00 0.94 0.68 
ECG Ep) 4.82 9.64 49.99 2067.00 0.94 1.27 
EC, 指数 3.05 6.27 48.61 1383.00 0.93 0.52 
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各 土 层 表 观 电导 率 的 基 台 值 介 于 6.27~10.79 之 间 ， 
表明 各 土 层 的 表 观 电导 率 均 存 在 着 显著 的 空间 变 
异 。 块 金 值 与 基 台 值 之 比 表 示 随 机 因素 引起 的 空 
间 蜡 质 性 占 总 空间 异 质 性 的 大 小 ,各 土 层 表 观 电导 
率 的 块 金 值 与 基 台 值 之 比 均 接近 于 50% ,表明 随机 
性 因素 和 结构 性 因素 共同 影响 研究 区 土壤 表 观 电 
导 率 的 空间 异 质 性 。 变 程 表示 空间 自 相 关 的 最 大 
距离 ,各 土 层 表 观 电导 率 的 变 程 介 于 1383~2579 m 
之 间 ,并 且 采 集 数据 时 ,各 样 点 距离 控制 在 300 m zz 
内 , 均 小 于 变 程 ,表明 合理 的 采样 距离 不 仅 可 以 提 
高 插值 的 精度 ,还 可 以 准确 地 揭示 研究 区 土壤 盐 渍 
化 的 空间 异 质 性 。 因 此 可 根据 土壤 盐 渍 化 空间 异 
质 性 划分 管理 分 区 。 
3.2 表 观 电导 率 与 光谱 指数 和 地 形 数据 的 相关 性 
分 析 

为 探究 表层 土壤 表 观 电导 率 (ECio3;) 与 地 形 数 
据 (DEM) 及 经 不 同时 期 遥感 影像 提取 所 得 光谱 指 


数 之 间 的 相关 性 ,将 DEM 数据 及 所 提取 的 22 个 光 
谱 数 据 与 表层 土壤 表 观 电导 率 数 据 在 Excel 中 进行 
相关 系数 计算 ,结果 如 表 3 所 示 。 通 过 分 析 棉 花生 
长 旺盛 期 9 月 的 9 个 植被 指数 可 知 ,DVI 和 GDVI 与 
土壤 表 观 电导 率 未 达到 显著 相关 ,而 RVI、.GRVI 和 
NDVI 等 7 个 植被 指数 与 表 观 电导 率 达 到 极 显著 相 
关 。 对 棉花 收获 后 冬 灌 前 11 月 裸 土 的 13 个 土壤 盐 
分 指数 分 析 可 知 ,NDSI、S3 和 S5 等 8 个 土壤 盐分 指 
数 与 土壤 表 观 电导 率 达 到 显著 相关 ,其 中 NDSI、S5 
和 SI-T 与 土壤 表 观 电导 率 达 到 极 显 著 正 相 关 。 此 
外 ,由 表 3 可 知 ,植被 指数 中 除 NG 外 其 余 均 与 表 观 
电导 率 呈 负 相 关 ,而 土壤 盐分 指数 却 呈 相反 结 

除 S6 外 其 余 均 与 表 观 电导 率 呈 正 相 关 ; 这 与 彭 杰 ” 
对 多 光谱 波段 和 光谱 指数 与 电导 率 的 相关 性 关系 
研究 结论 相 一 致 。 从 地 形 数据 分 析 可 知 ,DEM 数据 
与 表 观 电导 率 未 达到 显著 性 水 平 , 这 主要 原因 包括 
两 方面 ,一 是 研究 区 地 势 平 坦 ,削弱 了 地 形 条 件 对 


X3 表 观 电导 率 与 光谱 指数 和 地 形 数据 的 相关 性 分 析 


Tab.3 Correlation analysis table of apparent electrical conductivity with spectral index and DEM data 


光谱 指数 及 地 形 数据 相关 系数 计算 公式 

植被 指数 比值 植被 指数 (RVI) -0.11" NIR/R 
差 值 植被 指数 (DVI) -0.02 NIR-R 
绿 度 差 值 植 被 指数 (GDVI) -0.02 NIR-G 
绿 度 比 值 植被 指数 (CRVI) -0.12" NIR/G 
归 一 化 植被 指数 (NDVI) -0.09" (NIR-R)/(NIR+R) 
增强 植被 指数 (EVI) -0.11" 2.5x(NIR-R )/(NIR+6XR-7.5XB+1 ) 
近 红 外 归 一 化 (NNIR) -0.09” NIR/(NIR+R+G ) 
红 光 归 一 化 指数 (NR ) -0.09" R/(NIR+R+G) 
绿 光 归 一 化 指数 (NG ) 0.09” G/(NIR+R+G) 

土壤 盐分 指数 归 一 化 差 值 盐分 指数 (NDSI) 0.11" (R-NIR )/(R+NIR) 

盐 涡 化 指数 1(S1) 0.02 B/R 
盐 涡 化 指数 2(S2) 0.03 (B-R)/(B+R ) 
盐 涡 化 指数 3(S3 ) 0.07 (GXR)/B 
盐 渍 化 指数 5(S5 ) 0.10" (BxR)/G 
盐 溃 化 指数 6(S6 ) -0.03 (RXNIR )/G 
盐 渍 化 指数 (SI-T) 0.10™ (R/NIR)x100 
盐分 指数 (SI) 0.08" SQRT(BXR) 
盐分 指数 1(SI1) 0.07 SQRT(GXR) 
盐分 指数 2(SI2) 0.04 SQRT(G'2-R^24 NIR?2) 
盐分 指数 3(SI3 ) 0.07 SQRT(G2-R^2) 
强度 指数 1(Int1) 0.07 (G+R)/2 
强度 指数 2(Int2 ) 0.05 (G+R+NIR)/2 

地 形 数据 地 形 数据 CDEM ) -0.01 - 


波段 。 


注 :*** 表 示 通 过 0.01 .0.05 水 平 显 著 性 检验 。R、G、B、NIR 均 为 Lantsat8 波 段 反射 率 ,R 为 红 光 波段 ,G 为 绿 波段 ,B 为 蓝 波段 ,NIR 为 近 红 外 
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盐分 空间 分 布 的 文 配 作用 ,二 是 灌溉 .土壤 改良 等 
人 类 管理 活动 改变 了 地 形 条 件 控制 下 盐分 空间 分 
布 的 原始 格局 。 本 研究 在 相关 性 分 析 基 础 上 筛选 
出 相关 性 达到 极 显 著 , 且 相关 系数 rz0.10 的 植被 指 
数 和 土壤 盐分 指数 ,将 其 作为 精准 管理 分 区 的 辅助 
变量 ,可 有 效 提 高 分 区 内 同 质 性 和 分 区 间 异 质 性 。 


3.3 基于 表层 表 观 电导 率 与 不 同 变量 下 的 精准 管 
理 分 区 
基于 各 土 层 表 观 电 时 率 空间 结构 的 分 析 ,为 了 


更 加 直观 地 揭示 研究 区 土壤 盐分 的 空间 异 质 性 特 
征 , 使 用 普通 克 里 格 插值 法 绘制 出 各 土 层 的 表 观 电 
导 率 分 布 图 。 因 各 变量 间 单 位 不 同 , 需 将 筛选 出 的 
植被 指数 (RVI、GRVI、EVI) .土壤 盐分 指数 (NDSI、 
S5 .SI-T) 及 各 土 层 土壤 表 观 电导 率 插值 图 进行 标准 
化 处 理 ,将 变量 控制 在 0~1 之 间 。 人 然后 将 各 土 层 表 
观 电 导 率 插值 图 重 采 样 为 30 mx30 m 的 像 元 ,保持 
与 其 他 变量 像 元 大 小 一 致 。 在 棉花 收获 后 ,土壤 中 
大 部 分 的 盐分 在 生育 期 内 会 通过 蒸发 作用 将 深层 
盐分 随 水 运 移 到 表层 土壤 ,因此 本 文选 取 表 层 表 观 
电导 率 (ECosaxs) 为 主 变量 ,其 他 变量 为 辅助 变量 进 
行 管理 分 区 。 同 时 ,通过 eCognition 软件 提取 出 农 
田 边界 ,进而 进行 对 农田 精准 管理 分 区 。 

通过 多 次 多 尺度 分 割 试验 确定 最 为 合适 的 形 
状 因子 为 0.3, 紧 致 度 因 子 为 0.5。 在 控制 形状 因子 
与 紧 致 度 因 子 不 变 的 条 件 下 ,通过 最 优 分 区 评价 指 
数 曲 线 (图 2) 可 知 ,基于 ECwsws 的 最 优 分 割 尺 度 为 
19, 分 区 数 为 19; 基 于 EC 和 深层 表 观 电导 率 的 最 
优 分 割 尺 度 为 19, 分 区 数 为 16; 基 于 ECoas 和 植被 
指数 的 最 优 分 割 尺 度 为 14, 分 区 数 为 17; 基 于 
ECoss 和 土壤 盐分 指数 的 最 优 分 割 尺 度 为 14, 分 区 
数 为 13。 不 同 变量 下 的 管理 分 区 如 图 3 所 示 。 由 图 


(b) Z& T EC, 3754] 深层 


(a) SEF ECi s HUS E 表 观 电导 率 分 割 尺 度 
ior, 一 分 区 数 180 | 一 分 区 数 | 
140 — 评价 指数 pq 2.2 160 Hy 一 评价 指数 人 | 
120 20 x a 

x 100 ul R 100 
i 80 Ls z M gol 
R 60 1.6 E R 60 L 
40 40 上 
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5 1015202530 5 1015202530 


分 割 太 度 分 割 尺度 


ChinaXiv 合 作 期 刊 


2 期 白 建 铎 等 :基于 多 源 数据 的 干旱 区 盐 渍 化 农田 精准 管理 分 区 研究 651 


3 可 知 ,不 同 变量 下 的 分 区 结果 , 均 可 将 不 同 程度 电 
导 率 区 域 进 行 分 区 ,表明 基于 多 源 数 据 的 面向 对 象 
多 尺度 分 制 的 管理 分 区 具有 可 行 性 。 
3.4 不 同 输入 量 下 管理 分 区 对 比 与 评价 
由 图 3 对 比 基 于 不 同 数据 下 的 4 种 最 优 管理 分 

区 发 现 , 由 于 辅助 变量 的 差异 ,分 区 结果 也 存在 明 
显 差异 。 从 农田 的 耕作 角度 出 发 ,基于 ECioax 的 分 

区 结果 ,存在 相对 罕 而 长 的 分 区 ,这 种 分 区 与 实际 
机 械 化 工作 不 相符 ; 当 引 入 辅助 变量 后 ,分 区 结 
在 每 个 分 区 中 最 少 有 2-~3 个 田 块 ,与 实际 机 械 化 工 
作 条 件 相 符合 。 从 分 割 效果 出 发 ,基于 ECioas 的 分 
市 尺度 与 基于 ECiws 和 深层 表 观 电导 来 的 分 割 尺度 
相同 ,但 分 区 个 数 有 差异 ,这 可 能 是 由 深层 土壤 表 
观 电导 率 与 表层 表 观 电导 率 分 布 格局 不 同 引起 的 ; 
基于 EECoas 和 植被 指数 的 分 割 尺度 与 基于 ECias 和 
土壤 指数 的 分 割 尺度 相同 ,但 分 区 个 数 不 同 。 上 述 
对 比 只 是 对 不 同 分 区 结果 进行 了 定性 对 比 ,而 未 进 
行 定 量 评价 。 分 区 结果 的 定量 评价 , 需 本 着 分 区 内 
同 质 性 最 大 , 且 分 区 间 蜡 质 性 最 大 原则 展开 评价 ， 
即 分 区 内 平均 变异 系数 越 小 ,分 区 间 莫 兰 指数 数值 
越 接近 于 零 ,代表 分 区 效果 越 好 。 由 表 4 可 知 ,不 同 
数据 下 各 分 区 内 平均 变异 系数 在 13.60%~18.65% 之 
间 , 根 据 刘 宁 3 对 不 同 土地 利用 方式 下 黄河 三 角 洲 
土壤 特性 空间 变异 研究 可 知 ,变异 系数 在 0~20% 之 
间 属 于 弱 变 异 ,在 20%~50% 属 于 中 等 变异 , 5096 
属于 强 变 异 ; 说 明 各 分 区 内 均 属 于 弱 变 异 ; 且 较 整 
个 研究 区 的 变异 系数 下 降 了 约 60% ,说 明 分 区 内 同 
质 性 较 高 。 不 同 数据 下 各 分 区 间 莫 兰 指 数 均 接近 
于 0, 表 示 各 分 区 间 相 关 性 较 低 ,说明 各 分 区 间 异 质 
性 较 高 。 综 上 从 农田 耕作 角度 .分割 效 果 及 分 区 评 
价 原 则 出 发 ,基于 ECwws 和 深层 表 观 电导 率 的 管理 


(c) 基于 BCso37s 和 和 (d) 基于 ECnuoaxs 和 土壤 


植被 指数 分 割 尺 度 盐分 指数 分 割 尺 度 
] [ 一 分 区 数 
pvc 评价 指数 , -| 1.8 
| 4 L6 
1.4 = 
112 
1.0 = 
+ 0.8 
| 0.6 
0.4 -一 一 一 一 一 0.4 
4 8 12 16 20 4 8 12 16 20 
分 割 尺 度 分 割 尺度 


图 2 精准 管理 分 区 评价 指数 


Fig.2 Segmentation evaluation index 
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(c) 基于 BCsosws 和 植被 指数 的 分 区 结果 
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(b) 基于 ECio37s 和 深层 表 观 电导 率 的 分 区 结果 
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图 3 不 同 变量 下 的 最 优 管理 分 区 


Fig. 3 Optimal management zones under different variables 


表 4 基于 不 同 数 据 集 下 的 分 区 评价 


Tab.4 Partition evaluation based on different dataset 


分 区 信息 分 区 内 变异 系数 /% ”分 区 间 莫 兰 指数 
BCioss 13.60 0.1281 
ECnoss+ 深 层 表 观 电导 率 15.84 0.0075 
EC ss) UE BL 18.65 0.0184 
ECiosst THESE A) HR BL 16.30 0.0168 


分 区 为 最 佳 ,因为 该 分 区 结果 既 符 合 当 地 农田 的 管 
理 和 机 械 化 作业 要 求 ,又 可 以 代表 深层 盐分 的 分 
区 , 且 分 区 间 异 质 性 也 最 强 。 基 于 EC 和 土壤 盐 
分 指数 的 管理 分 区 次 之 。 基 于 ECiosx 的 管理 分 区 最 
差 , 是 由 于 其 分 区 结果 存在 着 罕 而 长 的 分 区 ,与 实 
际 耕 作 和 机 械 化 作业 不 相符 。 


4 讨论 


新 疆 南 疆 地 区 是 我 国 棉花 主要 生产 基地 ,但 该 
地 区 盐 省 化 造成 的 生态 环境 恶化 ,作物 生产 力 下 


降 .土壤 肥力 降低 问题 ,是 制约 该 地 区 作物 高 产 的 
主要 障碍 因子 ,因此 对 该 地 区 盐 渍 化 分 区 具有 重要 
意义 。 本 研究 选取 了 阿拉 尔 暴 区 棉花 收获 后 冬 灌 
前 的 农田 为 研究 对 象 ,此 时 期 盐分 聚集 在 土壤 表 
层 ” , 且 研 究 区 土壤 质地 以 砂 壤 土 为 主 , 保 水 保 肥 
能 力 差 ,土壤 中 水 分 ,养分 不 稳定 , 受 环境 影响 较 
大 ,因此 在 研究 区 进行 盐 污 化 分 区 比 进行 水 分 , 养 
分 分 区 更 有 意义 。 在 研究 方法 上 ,基于 面向 对 象 的 
多 尺度 分 割 算法 是 目前 进行 农田 管理 分 区 人 研究 中 
最 常用 的 方法 之 一 ,其 分 区 数 是 由 最 优 分 割 尺 度 所 
决定 ,不仅 消除 了 人 为 因素 的 影响 ,还 能 够 精准 地 
划分 管理 分 区 ,但 最 优 分 割 尺度 同时 受 形 状 因子 、 
紧 致 度 因 子 .影像 空间 分 辨 率 、 人 研究 区 面积 大 小 等 
多 个 因素 影响 。 本 文通 过 与 刘 焕 军 等 利用 遥感 
影像 对 黑土 区 田 块 尺度 的 精准 管理 分 区 对 比 发 现 ， 
影像 分 辨 率 不 同 其 最 优 分 割 尺度 也 不 同 , 且 最 优 分 
割 尺 度 随 着 空间 分 辨 率 的 增加 而 减 小 , 即 空间 分 辨 
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成 高 ,分 割 结 果 可 能 较为 破碎 ,其 分 区 结果 可 能 
适合 新 疆 这 种 高 度 机械 化 的 农事 作业 的 开展 。 
同时 最 优 分 割 尺度 与 研究 区 面积 大 小 也 具有 一 定 
关系 ,根据 黄 万 里 ”对 不 同 面积 的 森林 对 象 进行 多 
尺度 分 割 发 现 随 着 研究 区 面积 的 增 大 ,最 优 分 割 尺 
度 也 随 之 变 大 。 因 此 在 使 用 面向 对 象 多 尺度 分 割 
时 不 仅 要 考虑 人 研究 区 所 占 面积 大 小 ,更 要 选择 合适 
的 分 辨 率 ,从 而 确定 最 优 分 割 尺度 。 经 对 比 发 现 多 
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(3) 分 区 后 结果 表明 ,各 分 区 结果 的 变异 系数 
均 小 于 19% , 且 较 整个 研究 区 的 变异 系数 下 降 约 
60% ,说 明 面 向 对 象 多 尺度 分 割 算法 可 有 效 地 划分 
出 均 质 性 区 域 ,从 而 指导 农田 的 差异 化 管理 。 

(4) 对 比分 区 结果 发 现 ,基于 单一 数据 的 管理 
分 区 其 分 区 间 异 质 性 较 低 ,上 且 部 分 分 区 不 符合 耕作 
单元 ,而 基于 多 源 数据 的 管理 分 区 ,不 仅 符合 耕作 
单元 , 且 分 区 间 异 质 性 较 单一 数据 源 的 管理 分 区 均 
有 所 提高 。 

(5) 从 农田 的 耕作 角度 、 分 割 效果 及 分 区 评价 


源 数 据 较 单 一 数据 可 较 好 地 划分 管理 分 区 , 且 在 本 
质 上 更 具有 说 服 力 ”。 本 文 基于 辅助 变量 的 多 源 
数据 管理 分 区 消除 了 过 于 破碎 化 的 分 区 , 且 分 区 间 
的 异 质 性 也 有 明显 提高 ,但 可 能 存在 数据 元 余 和 计 
算 量 增加 等 问题 ,因此 在 后 期 研究 中 需 改 善 此 类 问 
题 。 在 分 区 管理 方面 , 接 下 来 将 根据 盐 渍 化 程度 计 
算出 各 个 分 区 压 盐 所 需要 的 灌溉 量 ,通过 测量 各 个 
分 区 的 土壤 自然 含水 量 .土壤 容重 ,田间 持 水 量 等 
基础 参数 ,采用 列 果 斯 塔 也 夫 经 验 公式 ,进行 各 个 
分 区 灌溉 量 的 计算 ” ,从 而 达到 精准 灌溉 的 目标 。 


5 结论 


本 研究 运用 以 电磁 感应 原理 为 基础 的 EM38- 
MK2 大 地 电导 率 仪 , 获 取 棉 花 收 获 后 冬 灌 前 表 观 电 
导 率 数据 。 用 半 变 异 函 数 进行 分 析 ,揭示 土壤 盐 渍 
化 的 空间 异 质 性 。 利 用 空间 分 辩 率 为 30 m 的 Land- 
sat8 影像 不 同 波段 信息 ,计算 出 22 个 光谱 指数 。 将 
光谱 指数 和 地 形 数据 分 别 与 表 观 电导 率 (ECoss) 展 
开 相 关 性 分 析 ,筛选 得 出 3 个 植被 指数 和 3 个 土壤 
盐分 指数 ,并 将 其 作为 精准 管理 分 区 的 辅助 变量 ， 
最 后 进行 多 源 数据 的 土壤 盐 渍 化 管理 分 区 研究 ,所 
得 结论 如 下 : 

(1) 在 半 变 异 函 数 分 析 中 ,各 土 层 表 观 电导 率 
的 基 台 值 与 块 金 值 之 比 均 在 50% 左 右 ,说 明 研 究 区 
土壤 盐 渍 化 存在 空间 变异 ,上 且 是 由 随机 性 因素 和 结 
构 性 因素 共同 导致 的 变异 ,因而 可 以 根据 研究 区 土 
壤 盐 渍 化 程度 进行 分 区 。 

(2) 经 相关 分 析 和 指数 筛选 后 ,9 月 植被 指数 
(RVI, GRVI, EVI) #il 11 月 裸 地 的 土壤 盐分 指数 
(NDSI.S5 SI-T) 均 与 ECoss 达 到 极 显著 水 平 ,因此 
将 其 作为 土壤 盐 渍 化 管理 分 区 的 辅助 变量 。 


原则 出 发 可 得 ,基于 ECwws 和 深层 表 观 电导 率 的 管 
理 分 区 为 最 佳 ,该 分 区 结果 既 符 合 当地 农田 管理 又 
符合 机 械 化 作业 ;而 基于 ECwws 和 土壤 盐分 指数 的 
管理 分 区 次 之 ,效果 最 差 的 为 仅 基 于 ECwoss 单 一 数 
据 源 的 管理 分 区 。 
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Precise management zoning in arid soil croplands based on multi-source data 
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Abstract: According to the spatial heterogeneity of soil salinization, the precise management precise zoning of 
oasis farmland in the arid area of Southern Xinjiang is of great significance for the adjustment of areas, based on 
spatial heterogeneity in soil salinization, is important for determining agricultural planting structure structures 
and fine-scale management. This study selected the topical oasis farmland in arid zone as the study object and 
used electromagnetic induction data, topographic date, and satellite remote sensing data as for oasis farmland in 
the data sources. We analyzed arid zone of Southern Xinjiang to analyze the spatial heterogeneity of soil 
salinization by geostatistical. Geostatistical methods and screened the were used to evaluate vegetation index and 
salt index in indices for different periods with correlation analysis. The surface apparent electrical conductivity 
(ECios;s) of the farmland was used as the main variable, and the deep apparent electrical conductivity (ECi;;, 
ECvo7s, ECyis), vegetation index (RVI, GRVI, EVI), and soil salinity index (NDSI, S5, SI- T) were used as the 
auxiliary variables. The study area was partitioned by using an object- oriented multi- scale segmentation 
algorithm, and the zoning results were evaluated by the mean Coefficient of variation (Cv) and Moran' s I. The 
results showed that there was obvious spatial heterogeneity in the apparent electrical conductivity of each soil 
layers in the study area, and the auxiliary variables were significantly correlated with the main variable. The 
average coefficient of variation of each partition was reduced by 6096 compared with that of the whole study 
area, and the. The interval heterogeneity based on multi-source data zoning was higher than that based on single- 
source data zoning. From the perspective of farmland cultivation, segmentation effects, and zoning evaluation 
principles, the, management zoning effect that integrated the information of surface and deep soil salinization was 
the best. And the zoning results were not only consistent with both local farmland management but also met the 
mechanized operation requirements. The findings can provide certain a technical and methodological reference 
for precise management zoning of oasis farmland in arid areas of Southern Xinjiang. 
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